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NOTATION 


A airfoil enclosed area divided by c 2 

a^ participation coefficient 

c chord length 

drag coefficient 
C lift coefficient 

li 

pitching-moment coefficient 
0^ pressure coefficient 

F(X) design objective function to be minimized or maximized 

f arbitrary function of X 

G.(X) constraint function 

J 

[H] Hession matrix containing second partial derivatives 

n number of design variables 

S q search direction in optimization 

t/c thickness-to-chord ratio 

X vector containing the design variables 

£ 

X_^ lower bound on design variable i 

X^ U upper bound on design variable i 

Y array of airfoil coordinates 

Y 1 array of coordinates defining a shape function 
a airfoil angle of attack 

a* move parameter in optimization 

V gradient operator 

A difference operator 



Subscripts : 


i variable number 

Is lower surface 

max maximum 

min minimum 

us upper surface 

Superscripts : 

0 nominal design 

q iteration number 

k design number 



APPROXIMATION CONCEPTS FOR NUMERICAL AIRFOIL OPTIMIZATION 


Garret N. Vanderplaats 
Ames Research Center 


SUMMARY 


An efficient algorithm for airfoil optimization is presented. The algo- 
rithm utilizes approximation concepts to reduce the number of aerodynamic 
analyses required to reach the optimum design. Examples are presented and 
compared with previous results. Optimization efficiency improvements of more 
than a factor of 2 are demonstrated. Much greater improvements in efficiency 
are demonstrated when analysis data obtained in previous designs are utilized 
The method is a general optimization procedure and is not limited to this 
application. The method is intended for application to a wide range of engi- 
neering design problems. 


INTRODUCTION 


Numerical optimization techniques have been shown to provide a versatile 
tool for airfoil design. The usual approach has been to couple existing aero- 
dynamic analysis codes with an optimization code to achieve the design capa- 
bility. The primary effort has been directed toward application of these 
techniques to a wide variety of design problems while using increasingly 
sophisticated and time-consuming aerodynamic analysis programs. 

The cost of this automated design, whereby a very time-consuming analysis 
program is executed repetitively (perhaps several hundred times) , is neces- 
sarily an important consideration when judging the practicality of these 
techniques. Perhaps the simplest means of estimating cost is by the number of 
times the aerodynamic analysis program is executed during a design study. 

That is, for a given aerodynamics program, the cost of optimization is a 
direct function of the number of times the program is executed on the com- 
puter. Therefore, any improvement in optimization efficiency is directly 
measurable in design cost savings. 

Very little effort has been directed toward improving the efficiency of 
the automated design process as applied to aerodynamic design. The principal 
improvement to date has been in the method of defining the airfoil. In refer- 
ences 1 and 2, polynomials were used to define the airfoil shape, with the 
coefficients of the polynomial being the design variables. In references 3 
and 4, and in subsequent work, these polynomials were replaced by more general 
analytical or numerically defined shape functions. The result was an effi- 
ciency improvement of more than a factor of 2, together with improved airfoil 



definition (ref. 3). However, efficiency improvements are still needed if 
numerical airfoil optimization is to become an economically feasible design 
approach when using sophisticated aerodynamic analysis codes. 

The purpose here is to present a technique that improves the design effi- 
ciency by another factor of 2 or more. The basic approach is to develop 
approximations to the design problem using a minimal amount of information. 

The approximating functions are used in the optimization and the resulting 
design is analyzed precisely. This analysis information is added to the 
available data and the process is repeated until convergence to the optimum 
is achieved. (The idea of using approximation concepts in aerodynamic opti- 
mization originated from the observed success of similar techniques used by 
Schmit and Miura (ref. 5) in the field of structural optimization.) 

To provide a background for the method, the basic concepts of numerical 
optimization are first presented. This is followed by a description of the 
present method of coupling an aerodynamic analysis program to an optimization 
program for automated design. The concept of optimization by sequential 
approximations is then presented, followed by a more precise mathematical 
formulation of the method and a summary of the design algorithm. Examples 
demonstrate the efficiency and reliability of the method and finally, some 
of the implications of this method for future development are discussed. 


OPTIMIZATION CONCEPTS 


Assume the airfoil is defined by the relationship 

Y = a, Y 1 + a 2 Y 2 + . . . + a Y n 
1 z n 


( 1 ) 


where Y is a vector containing the upper and lower coordinates of the air- 
foil and yi are shape functions that may themselves define airfoils. The 
coefficients a^, a 2 . . . , a n are referred to as participation coefficients. 
Now assume it is desired to find the airfoil that minimizes the drag coeffi- 
cient Cp with constraints on lift coefficient C L , thickness-to-chord ratio, 
t/c, etc., at a specified Mach number and angle of attack. The participation 
coefficients aj - a R are the design variables, and will be changed during 
the optimization process. The n-dimensional space spanned by the design 
variables is called the design space. 

The optimization problem can be stated mathematically as: 


Minimize 


subject to 


C 


D 


C L > 


"^mi 


m 


( 2 ) 

(3) 


2 



(4) 


(t/c) > (t/c) 


min 


where C^, C^, and t/c are functions of a 1? a 2 , . . . , a n . This can be 
generalized to be: 

Minimize 


F(X) 


subject to 


G.(X) < 0 j = 1, m 

x. £ < X. < X. u i = 1, n 
111 


(5) 

( 6 ) 
(7) 


where X is a vector containing the design variables, a^, i = l,n. There are 
a total of m constraints. The lift constraint of equation (3) is written in 
the form of equation (6) as 


1 



( 8 ) 


Similarly, from equation (4) 


1 


(t/c) 


(t/ c) 


< 0 


mm 


(9) 


it U 

The parameters X^ and Xj_ of equation (7) are referred to as side con- 
straints that limit the region of search for the optimum. Although side 
constraints could be included in equation (6) , they are usually treated sepa- 
rately for convenience and efficiency. 


The optimization problem of equations (5)-(7) is quite general. If it 
is desired to maximize with a constraint on Cq, -C^ is minimized. 

Also, the constraint set of equation (6) is not limited to constraints at the 
design flight condition. With this formulation, the airfoil can be designed 
at a given flight condition with constraints at other flight conditions so 
long as the appropriate information is calculated during the aerodynamic 
analysis. If the inequality conditions of equations (6) and (7) are satis- 
fied, the design is said to be feasible. If any of these conditions are 
violated, the design is called infeasible. 


The optimization process typically proceeds in an iterative fashion as: 

x q = X q_1 + a* s q (10) 

An initial design, X° , must be provided which may or may not define a feasible 
design. The superscript q is the iteration number. Vector is the 
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search_direction and a* is a scaler determining the move distance in direc- 
tion sq. The notation a is used for consistency with mathematical pro- 
gramming literature and should not be confused with the airfoil angle of 
attack. 


If gradient methods are used, the optimization process consists of two 
steps. The first is determination of a move direction that will improve 

the design without violating any constraints; the second is calculation of 
a* such that the objective is reduced as much as possible in this direction. 


This may be understood by considering a two-variable design problem where 
C D is minimized subject to constraints on Cl and t/c. A hypothetical 
problem is shown in figure 1 which shows contours of constant drag and shows 
the Cl and t/c constraint boundaries. Assume an initial design is given at 
point A, with no active or violated cons traints . _ Using gradient methods, the 
process begins by perturbing each component of X to determine its effect on 
the objective. That is, the gradient of Cp is calculated by finite differ- 
ence using a single forward step, and the gradient vector is constructed as 


3c d ] 


( ac d ) 

dX 1 


&X 1 

3c d 


ac d 

3X 2 


ax 2 


► ~ < 


9C d 


ac d 

ax 

n J 


AX 

n / 


(ID 


It is obvious that the greatest improvement of the objective function is 
achieved by moving in the negative gradient, or steepest descent direction, 
so that S = -VCq. Knowing the search direction, S 1 , the scalar a* that 
will minimize C^ in this direction must be found. This is a one-variable 
minimization problem. Several somewhat arbitrary values of a are defined 
and the airfoil is analyzed at each point, X = X° 4- aS 1 . A polynomial is 
usually fit to these points and a more precise a = a* is calculated at point 
B in figure 1, ending the first optimization iteration. The second iteration 
begins by again perturbing the design variables to obtain the gradient of the 
objective. Using the conjugate direction algorithm of Fletcher and Reeves 
(ref. 6) a new search direction, S 2 , is found which will again reduce the 
objective. A search is performed in this direction, leading to point C, 
completing iteration two. At C, the lift constraint is active (Gj - 0) and 
a direction is found that will reduce Cp without violating this constraint. 
The gradient of both the objective and active constraint are calculated and 
a new search direction, S 3 , is found using Zoutendijk’s method of feasible 
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directions (refs. 7,8). The process is repeated until a design at E is 
obtained where no direction can be found that will reduce the objective with- 
out violating the constraints and this design is called optimum. Logic is 
included in the algorithm so that if an initial design is defined at point F, 
the constraint violations are overcome to yield a final design at point E. 

The optimization procedure described above is essentially that used in 
the CONMIN program (ref. 9). In a typical design, about 10 iterations are 
required to achieve an optimum design. For each iteration, n aerodynamic 
analyses are used to calculate the required gradient information by finite 
difference. To determine a* requires an average of three analyses so that 
a total of 10(n + 3) = lOn + 30 aerodynamic analyses are required for opti- 
mization for a single flight condition. Although quite efficient from an 
optimization viewpoint, that many executions of a sophisticated aerodynamics 
program can be very expensive. Therefore, it is desirable to reduce the number 
of required analyses as much as possible. This improvement in optimization 
efficiency is the subject here. The technique will be developed by first 
reviewing the approach currently used for aerodynamic optimization. 


Previous Method of Aerodynamic Optimization 

At the present time, most airfoil optimization is performed by coupling 
the aerodynamics program to the optimization program as shown in figure 2. 

Each time the optimization program defines a new design, either for finite 
difference gradient computations or for determining a*, the aerodynamics 
program is called for a complete analysis. For the example of figure 1, a 
set of analyses is performed as indicated in figure 3. During optimization, 
at iteration j, very little information from previous iterations is used. 

At point B in figure 1, the vector S 1 is used to calculate S 2 so that, if 
no constraints are active, prior information is used. However, if one or 
more constraints are active or violated (the usual situation) , no prior 
information is used. 

It can be argued intuitively that all calculated information should 
be of value in guiding the optimization process. Furthermore, in a design 
study, numerous optimizations are usually performed. For example, one 
optimization may be done to minimize C D with constraints on and, later, 

another optimization done to minimize with constraints on and C^. 

It may be expected that, because many airfoils were analyzed during the first 
optimization, a second optimization at the same flight condition should 
utilize this available information. One way to do this would be to approxi- 
mate the required functions using available information. This would provide 
explicit functions which could now be optimized independent of the time 
consuming aerodynamic analysis program. Aerodynamic analysis is still used 
to improve the approximation, leading to a precise solution. 

The general procedure is outlined in the following section. 
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Optimization by Sequential Approximations 


Assume it is desired to approximate the aerodynamic parameters linearly 
in terms of the design variables, X. Then, for example. 


3C 

r ~ r 0 j 

S D 3X X 


ac 

(x i - x i 0) + nr 


(X 2 - X 2 °) + . 


3C. 


D 


9X 


n 


(X - X °) 
n n 


( 12 ) 


= C D ° + VC D • (X - X°) (13) 

where the superscript denotes the point about which the curve fit was done. 

Now at point A in figure 1, this information is obtained by n + 1 analyses; 
an initial function evaluation plus finite difference gradients. Similar 
information can be simultaneously obtained for Cl and t/c. This provides 
explicit, but approximate, expressions for the functions. The two-variable 
representation for this linearized problem is shown in figure 4, where the 
objective and constraints are now linear functions of the design variables. 

It is clear from figure 4 that the solution to the approximate optimization 
problem is unbounded so that this would not yield meaningful results. 

However, by limiting the design change to some reasonable bounds, shown by 
the rectangle, the optimization will have a solution. At that point, the 
problem could be linearized again and resolved. Note that at the end of this 
first optimization the airfoil will be analyzed precisely and the results 
compared with the approximate solution. Therefore, n 4- 2 analyses are now 
available to set up a new approximate problem. The problem could be linear- 
ized by either using only n + 1 of these analyses or using a least squares 
fit to all of the data. 


Note that, because this approximate problem is strictly linear, linear 
programming techniques could be used (Simplex method, ref. 10). The technique 
of repeatedly linearizing a nonlinear problem and solving with linear pro- 
gramming is known as "sequential linear programming" and has been used with 
success in structural design (ref. 11). 

An alternative approach to sequential linear programming, and the one 
used here, is to use the excess data to develop higher order approximations 
to the functions. Then, in this case, the extra analysis would be used to 
provide a second-order approximation with respect to Xj . This new approxi- 
mate problem is optimized, followed by a new precise analysis and the process 
repeated until the solution has converged. When a total of 

1 + n 4- n(n + l)/2 analyses are available, a full quadratic approximation is 
possible. Subsequent to this, new analyses are used in a weighted least 
squares fit, rather than obtaining ever higher order approximations. Only a 
second-order approximation is used because higher order approximations 
(1) would require excessive data, (2) tend to model noise in the data, and 
(3) have been found to be unnecessary. 

If a quadratic design problem were being considered, the quadratic 
approximation would be precise. If the problem at hand can be approximated 
closely by a quadratic function, this method can be expected to be 
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competitive. Assuming only a few analyses are required beyond that necessary 
for a quadratic approximation, this method is competitive for problems of 
fewer than 20 design variables, and is twice as efficient for problems of 
10 variables. More dramatic improvements are realized by using analysis data 
obtained in one optimization as data for subsequent optimizations. 

The concept of sequential approximations is shown in block diagram form 
in figure 5. Note, in comparison to figure 2, that the optimization program 
never directly works with the aerodynamic analysis, but only optimizes the 
approximating functions. Because the evaluation of the approximate functions 
is short and is explicit, the efficiency of the optimization code itself is of 
minor importance because the necessary function evaluations are quite rapid. 
Also, gradients of the approximating functions are easily calculated analy- 
tically, a feature that improves the efficiency of the approximate 
optimization. 

In the following section, the sequential approximation approach to 
optimization is outlined mathematically. 


Mathematical Formulation 


Consider the Taylor series expansion of any function: 


f = f° + Vf • (X - X°) + (X - X°) • H • (X - X°) + • • • (14) 

In equation (14) , X° is the point about which the expansion is being per- 
formed, f° is the corresponding function value, and Vf is the vector of 
first partial derivatives (gradient). The matrix of second partial deriva- 
tives (Hession matrix), H, is symmetric: 


9f 


9 2 f 9 2 f 

3X X 


9X^2 9X X 9X 2 

3f 


9 2 f 

3X 2 


9X 2 2 

• 

► H = « 


3f 

9X n 

k i 


SYMMETRIC 


3 2 f 

3x 1 sx n 


(15) 


3 2 f 
3X n 2 


/ 


In aerodynamic optimization, this information is not usually available analy- 
tically, but can be calculated numerically. The usual approach is to use 
finite difference gradients to provide a good approximation to Vf and H. 

This assumes a small finite difference step size to insure that the 
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approximation is good. However, it may be desirable to utilize previously 
calculated data to determine the components of Vf and H, recognizing that the 
accuracy of this approximation will not be as great as if small finite differ- 
ence steps were used. 

The terms up through second order of equation (14) may be written in more 
compact form as 

Af = Vf • AX + j AX T [H] AX (16) 


where 


AX = X - X° 


and 


Af = f - f° 


Equation (16) may be expanded as 


Af = Vf] AX-, + VfoAXo + • • • Vf AX 
11 z ^ n n 



AX X 2 + H 2 2 AX 2 2 + 


+ H AX 2 ) 
nn n 


+ H^AXxAXa + E 13 AX l AX 3 + • ■ • 4- H] n AXi AX n 


+ h 23 ax 2 ax 3 + • • 

Now assume a nominal design, X° , has 
numerous other designs, X 1 , . . . X^ 

f 1 , . . . , fk. Let 

AX 1 = X 1 - X° 

and 

Af 1 = f 1 - f° 

Then we can write k equations of the form (17) . The unknowns are 
Vfi . . . Vf n and Hj i , H 12 . . . H nn for a total of l = n + n(n 4- l)/2 
unknowns. Remembering that one analysis was required for the nominal design, 
X°, a total of i + 1 designs are required. Thus, if k > £ + 1, the 
unknowns can be determined. Equation (17) is linear in the unknown coeffi- 
cients. Writing this equation for each of l designs yields £ equations 
which can be solved directly. 

If more than £ designs are available, a weighted least squares fit is 
used. If less than £ designs are used, fewer coefficients are calculated. 


’ H n-l,n Ax n-1 AX n ( 17 ) 

been analyzed to yield f°. Also, 
have been analyzed to give 

i = 1, - • • » k 
i = 1 , . . . , k 
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In the extreme, if only the nominal design and one other design are available, 
only the first term in equation (17) is calculated. This approximation can 
be used to optimize with respect to variable only. The result of that 

optimization is then analyzed and used to calculate the first two terms in 
equation (17). Optimization can then be performed with respect to and 

X 2 , and so on. In this fashion, all data are used to guide the optimization, 
and the design is continually improving. Figure 6 depicts the sequence of 
designs that may be analyzed precisely in the two variable example. 


Design Algorithm 

Given the capability of developing the approximate Taylor series expan- 
sion of the various aerodynamic and geometric functions, it is incorporated 
into an optimization algorithm as follows: 

1. Given k initial designs, k > 2 

2. Create the Taylor series expansion about the current "best" design 

3. Number of design variables, NDV = MIN(k,n) 

4. Set limits on the design variables, say X^ = 0.8*X° and 
X u = l.2*X° 

5. Optimize the approximating functions 

6. Analyze the proposed optimum 

7. Add results to data set; set k = k + 1 

8. If k < n + 1 go to step 2 

9. Check convergence 

10. If satisfied, print final results; otherwise, go to step 2 

A FORTRAN computer code was written for this technique and a block 
diagram of the major operations is shown in figure 5. 

The CONMIN program (ref. 9) was used for the optimization capability. 

In the following section, design examples are presented to demonstrate 
the efficiency of the method. 


DESIGN EXAMPLES 


Examples are presented here to identify the generality and efficiency 
of approximation concepts as applied to airfoil optimization. Four existing 
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airfoils are used as the design basis vectors in equation (1). These are the 
NACA 2412, NACA 64 i -412, NACA 652-415, and the NACA 64 2 -A215 airfoils. The 
coordinates are defined at 50 points along the upper and lower surfaces. The 
coordinates are approximate, obtained from curve fits of the existing airfoil 
data (refs. 12, 13) and no attempt was made to precisely match the data given 
in the references. Two additional basis vectors were used to impose the 
boundary conditions at the trailing edge of the airfoil: Y us = X/C, Y^ s = 0 
and Y us = 0, Y^ s = -(X/C). The shapes defined by these six basis vectors are 
shown in figure 7. These basis vectors are the same as those used in 
reference 3. For consistency, the same aerodynamic analysis code (ref. 14) 
was also used. Three of the design examples of reference 3 are solved here, 
two of which were also presented in reference 2. 


Example 1: Lift Maximization, M = 0.1, a=6° 

Figure 8 shows the results of optimization of an airfoil for maximum 
lift. The design constraints are listed on the figure and are the same as 
Case 2 of reference 3. Additional numerical results are given in tables 1 
and 2. This optimization required 19 aerodynamic analyses; 44 were required 
previously. Although it may be, argued that this airfoil is impractical, it 
must be remembered that this airfoil mathematically satisfies the design 
constraints. Also, the lift coefficient obtained here, Cl = 1.144, is better 
than the one obtained before, Cl = 1.108. The fact that this airfoil was not 
obtained using the method of reference 3 suggests that the present method is 
numerically better conditioned for optimization. Furthermore, these results 
were obtained using less than half the number of aerodynamic analyses used 
previously. 

The quality of the approximation to the lift coefficient may be judged 
from figure 9. Because there are four independent design variables, the 
full second-order Taylor series expansion of the functions requires 15 anal- 
yses. It is intriguing to note that on the sixteenth analysis and beyond, 
the approximation for this case is quite precise. 


Example 2: Lift Maximization, M = 0.75, a = 0° 

Although optimization using sequential approximations works well for 
low-speed airfoils, it might be expected that the technique would not be 
adequate for high-speed applications where the nature of the flow field about 
the airfoil can be quite sensitive to small changes in the airfoil shape. 

To study this possibility, Case 3 of reference 3 was solved using the present 
method. Here, the lift coefficient was maximized subject to a constraint 
on wave drag. A value of Cl = 0.4211 was obtained after 27 analyses, 
compared to C L = 0.4188 obtained in 70 analyses previously. The results 
are shown as example 2A in figure 10. Using the present method, the optimi- 
zation continued to mathematically improve the airfoil. After 37 analyses, 
the airfoil of figure 11 (example 2B) was obtained. The optimization 
terminated after 48 analyses, yielding the airfoi] shown in figure 12 
(example 2C) . Note the significant changes in pressure distribution among 
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figures 10-13. The comparison of approximate and precise values of and Cq 

are shown in figure 13. The two values agree well at 27 analyses. After 37 
analyses the flow field is as shown in figure 11. The comparison between 
approximate and precise values is then poor until the convergence to the final 
optimum after 48 analyses. The disagreement between approximate and precise 
function values is understandable from figure 14 which shows the airfoil 
corresponding to the forty-third analysis. Note the reverse curvature of the 
upper surface near the leading edge and near the 60-percent chord. This 
results from the inability to properly model a supercritical airfoil using 
the NACA basis airfoils. The optimization was able to effectively utilize 
these data to redirect the optimization process, leading to the final 
converged solution. This final airfoil, shown in figure 12, has the same 
characteristics as the airfoil of figure 14, but to a lesser degree. In a 
practical design situation, it would be desirable at this point to add other 
basis vectors that represent supercritical airfoils, remembering that the 
48 analyses already obtained provide useful data for the expanded optimiz- 
ation. The results of this optimization are given in tabular form in tables 1 
and 2 . 


Example 3: Wave Drag Minimization, M = 0.75, a = 0° 

To demonstrate the efficiency of the present method when multiple 
optimizations at the same flight condition are performed, Case 4 of 
reference 3 was solved. The 48 analyses performed to solve example 2 here 
were used as initial data. In reference 3, the optimum airfoil from the 
previous design was used as a starting point for this design. In the present 
study, the twenty-seventh analysis (fig. 10) was used as the nominal design 
about which the first Taylor series expansion was performed. An optimum 
design of Cp - 0.0009 was obtained using only two additional analyses. The 
resulting airfoil is given as example 3A in figure 15 and in tables 1 and 2. 
This result compares to an optimum Cp - 0.0007 obtained previously using 
44 aerodynamic analyses. 

As an additional exercise, this design was repeated beginning with the 
forty-eighth analysis of example 2 as the initial nominal airfoil. An 
optimum Cp = 0.0003 was obtained using four additional analyses. This 
design (example 3B) is presented in figure 16 and in tables 1 and 2. As 
seen from the figures, examples 3A and 3B represent quite different airfoils, 
although the actual calculated wave drag is negligible in each case. 


DISCUSSION 


Approximation concepts as applied to aerodynamic design have been 
presented. The technique has been shown to be more versatile and efficient 
than earlier techniques. It is not limited to two-dimensional airfoils, or 
to a single flight condition, and it is not limited to the aerodynamic 
analysis code used to provide the examples. The technique is a general 
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automated design procedure that may be applied to a wide variety of engineer- 
ing design problems in addition to the one considered here. 

The results presented here are considered preliminary and, as experience 
is gained through application of the method, further refinements can be 
expected. 

Future effort will concentrate on development of the computer code, 
written as part of this study, into a generally available code for distri- 
bution, applicable to problems of broad engineering interest. Efforts in 
aerodynamic design will be directed toward extension to more general design 
situations. Of fundamental importance is development of data storage and 
retrieval systems so that the ever increasing body of available aerodynamic 
data can be easily utilized in design. Finally, effort will be directed 
toward the use of experimental data as a basis for design. The general goal 
is to develop a distributable computer program and data base that the user 
can apply to his particular design problem at extremely low cost. 


Ames Research Center 

National Aeronautics and Space Administration 

Moffett Field, California 94035, November 17, 1978 


12 



REFERENCES 


1. Hicks, R. M. ; Murman, E. M. ; and Vanderplaats , G. N. : An Assessment of 

Airfoil Design by Numerical Optimization. NASA TM X-3092, July, 1974. 

2. Vanderplaats, G. N. ; Hicks, R. M. ; and Murman, E. M. : Application of 

Numerical Optimization Techniques to Airfoil Design. NASA Conference 
on Aerodynamic Analyses Requiring Advanced Computers, Langley Research 
Center, Va. , NASA SP-347, Part II, 1975. 

3. Vanderplaats, G. N. ; and Hicks, R. M. : Numerical Airfoil Optimization 

Using a Reduced Number of Design Coordinates. NASA TM X-73,151, 1976. 

4. Hicks, R. M. ; and Vanderplaats, G. N. : Application of Numerical 

Optimization to the Design of Supercritical Airfoils Without Drag 
Creep. 1977 SAE National Business Aircraft Meeting and Exposition, 
Wichita, Kans., 1977. 

5. Schmit, L. A.; and Miura, H. : Approximation Concepts for Efficient 

Structural Synthesis. NASA CR-2552, 1976. 

6. Fletcher, R. ; and Reeves, C. M. : Function Minimization by Conjugate 

Gradients. British Computer Journal, vol. 7, no. 2, 1964, pp . 149-154. 

7. Zoutendijk, G. G. : Methods of Feasible Directions. Elsevier, 

Amsterdam, 1960. 

8. Vanderplaats, Garret N. ; and Moses, Fred: Structural Optimization by 

Methods of Feasible Directions. Journal of Computers and Structures, 
vol. 3, 1973, pp. 739-755. 

9. Vanderplaats, Garret N. : CONMIN — A FORTRAN Program for Constrained 

Function Minimization: User’s Manual, NASA TM X-62,282, 1973. 

10. Fox, Richard L. : Optimization Methods for Engineering Design. Addison- 

Wesley, Mass., 1971. 

11. Moses, F. : Optimum Structural Design Using Linear Programming. Proc . 

A.S.C.E., 90, ST6, 1964, pp. 89-104. 

12. Abbott, Ira H. ; von Doenhoff, Albert E. ; and Stivers, Louis S., Jr., 

Summary of Airfoil Data. NACA Rep. 824, 1945. 

13. Abbott, Ira H. ; and von Doenhoff, Albert E. : Theory of Wing Sections. 

Dover Publications, Inc., New York, 1959. 

14. Jameson, Anthony: Transonic Flow Calculations for Airfoils and Bodies 

of Revolution. Grumman Aerodynamics Report 390-71, 1971. 


13 



TABLE 1.- DESIGN INFORMATION 





Coefficients 


Example 

Comment 


.. - . 





a l 

a 2 

a 3 

a 4 

i 

Initial design 

1.0 

0 

0 

0 

1 

Optimum design 

4.9323 

-0.6653 

-0.4802 

-2.3018 

2,3 

Initial design 

1.0 

0 

0 

0 

2A 

After 27 analyses 

.4989 

.4004 

-.5997 

.6181 

2B 

After 37 analyses 

.0533 

.9713 

-1.3342 

1.2632 

2C 

Optimum design 

-.4518 

1.3518 

-1.8513 

1.8959 

2D 

Forty-third analysis 

-.4954 

1.2148 

-1.8393 

2.0302 

3 A 

Initial analysis no. 27 

.4205 

.1930 

-.6256 

.8698 

3B 

Initial analysis no. 48 

-.2137 

-.7223 

-1.3723 

1.7285 
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TABLE 2.- ANALYSIS INFORMATION 


Example 

C L 

C Dw 

C M 

A a 

t/c 

Camber 

Number of analyses 

Ref. 2 

Ref . 3 

New 

method 

1 init. 

0.9864 

0 

-0.0645 

0.0809 

0.1193 

0.0201 

— 

_ 

— 

1 opt. 

1.1437 

0 

-.7496 

.0750 

.1501 

.0398 

103 

44 

19 

2,3 init. 

.4824 

.0113 

-.1068 

.0821 

.1201 

.0201 

- 

- 

- 

2A 

.4211 

.0039 

-.1017 

.0750 

.1119 

.0101 

143 

42 

27 

2B 

.4702 

.0039 

-.1341 

.0750 

.1154 

.0108 

143 

70 

37 

2C 

.4716 

.0040 

-.1505 

.0755 

.1191 

.0117 

143 

70 

48 

2D 

.4024 

.0019 

-.1319 

.0759 

.1191 

.0124 

- 

- 

- 

3A 

.3057 

.0009 

-.0770 

.0751 

.1108 

.0077 

- 

44 

2 

3B 

.3009 

.0003 

i 

-.0945 

.0765 

.1168 

.0079 

- 

44 

4 


a Cross-sectional area divided by the chord squared. 
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Figure 6.— Sequence of precise analyses using approximation techniques. 
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BASIS SHAPE 4 



BASIS SHAPE 5 



(b) Shapes 4 through 6. 


Figure 7.— Concluded. 



CONSTRAINTS: |C n (x/c = 0.01)1 < 2.0 |C M |< 0.075 A >0.075 

p us v 

t/c <0.15 CAMBER <0.04 



AIRFOIL SHAPE 



Figure 8.— Example 1: Lift maximization, M = 0.1, a 



Figure 9.— Optimization history for example 1. 
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